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ABSTRACT 

The basic cosmological distances are linked by the Etherington cosmic dis- 
tance duahty relation, r]{z) = Dl{z){1 + z)~'^ /Da{z) = 1, where Dl and Da 
are, respectively, the luminosity and angular diameter distances. In order to test 
its validity, some authors have proposed phenomenological expressions for r]{z) 
thereby deforming the original Etherington's relation and comparing the result- 
ing expressions with the available and future cosmological data. The relevance of 
such studies is unquestionable since any violation of the cosmic distance duality 
relation could be the signal of new physics or non-negligible astrophysical effects 
in the usually assumed perfectly transparent Universe. 

In this letter, we show that under certain conditions such expressions can be 
derived from a more fundamental approach with the parameters appearing in the 
ri{z) expression defining the cosmic absorption parameter as recently discussed by 
Chen and Kantowski. Explicit examples involving four different parametrizations 
of the deformation function are given. Based on such an approach, it is also 
found that the latest Supernova data can also be explained in the framework of 
a pure cold dark matter model (Einstein-de Sitter). Two different scenarios with 
cosmic absorption are discussed. Only if the cosmic opacity is fully negligible, 
the description of an accelerating Universe powered by dark energy or some 
alternative gravity theory must be invoked. 

Subject headings: Supernovae:general, cosmic distances, duality relation 
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1. Introduction 

The cosmic distance duality relation ( CDDR ^ is a mathematical identity relating the 
luminosity distance Dl with the angular diameter distance Da by the expression: 

^{l + z)-'-l. (1) 

The validity of this constraint uniting the two basic distances in cosmology depends neither 
on the Einstein field equations nor on the nature of the matter-energy content. It only re- 
quires the phase space conservation of photons and that sources and observers are connected 
by null geodesies in a Riemannian spacctimc. Therefore, it remains valid for spatially homo- 
geneous and isotropic (anisotropic) cosmologies, as well as for inhomogeneous cosmological 
models (Etherington 1933, Basset and Kunz 2004). 

The above relation is usually taken for granted when relativistic models of the Universe 
are confronted to the existing cosmological observations. Despite that, the distance-duality 
relation is in principle testable by means of astronomical observations. The basic idea is 
to find cosmological sources whose intrinsic luminosities are known (standard candles) as 
well as their intrinsic sizes (standard rulers). After determining both Dl and Da at the 
same redshift it should be possible to test directly the Etherington result. Naturally, by 
cosmological sources with known D^ and Da we are not referring necessarily to the same 
class of objects. Under certain conditions, as recently discussed by Holanda et al. (2010), one 
may consider two different classes of objects as, for instance, supernovae and galaxy clusters 
for which D^ and Da are separately determined. Note also that ideally both quantities 
must be measured in such a way that any relationship coming from specific cosmological 
models are not used, that is, they must be determined by means of intrinsic astrophysically 
measured quantities. In practice, the vahdity of the CDDR has been tested by assuming a 
phenomenological deformed expression of the form (Holanda et al. 2010, Meng et al. 2011, 
Nair et al. 2011, Khedekar & Chakaborti 2011, Gongalvez et al. 2011): 

^il + z)-'^r,{z), (2) 

where ri{z) is the deformation function which quantifies a possible epoch-dependent departure 
from the standard photon conserving scenario [f] = 1). 

As it appears, a deformed CDDR can also be adopted to test the possibility of a new 
physics. In this line. Basset & Kunz (2004) used supernovae type la (SNc la) data as 
measurements of the luminosity distance and the estimated Da from FRIIb radio galaxies 
(Daly & Djorgovski 2003) and ultra compact radio sources (Gurvitz 1999; Lima & Alcaniz 
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2000, 2002; Santos & Lima 2008) in order to test possible new physics signatures based on 
the following expression: 

/. 7]{z) = (l + ^)^-^exp 

where E{z') = H{z')/Hq is the dimensionless Hubble parameter {Hq is the Hubble constant). 
Note that for arbitrary values of 6, the strict validity of the DD relation corresponds to 
(/3,7) = (1,0). By marginalizing on VIm, and Hubble parameters, they found a 2a 
violation caused by excess brightening of SNIa at 2; > 0.5. It was also argued that such an 
effect would be associated to lensing magnification bias. 

De Bernardis, Giusarma & Melchiorri (2006) also searched for deviations of the standard 
CDDR by using the angular diameter distances from galaxy clusters provided by the sample 
of Bonamente et al. (2006). By assuming 77(2;) = constant, they obtained a non violation 
of CDDR in the framework of the cosmic concordance ACDM model. Later on, Avgoustidis 
ct al. (2009, 2010) working in the context of a flat ACDM model also adopted an extended 
CDDR expressed as 

//. 7]{z) = {1 + zy. (4) 

The above deformation function is a particular case (/3 = 1 + e, 7 = 0) of the general 
expression adopted by Basset and Kunz (2004). In their analysis, the recent SNe la data 
as compiled by Kowalski et al. (2008) were combined with the latest measurements of the 
Hubble expansion at redshifts in the range < 2; < 2 (Stern et al. 2010), and the free 
parameter was constrained to be e = — 0.041q;q7 {2a). More recently, such a parametrization 
has also been adopted by Khedekar & Chakaborti (2011) in their studies connecting the 
Tolman test and CDDR using the redshifted 21 cm wavelength from disk galaxies. It was 
also argued that future data from the planned Square Kilometer Array (SKA) my provide 
the best test to detect any violation of the cosmic distance duahty relation. 

In a series of papers, Holanda, Lima & Ribeiro (2010, 2011, 2011a) also explored a 
different route to test the CDDR based on the following deformation functions: 

///. ri{z) = 1 + rjoz, IV. r]{z) = 1 + r]oz/{l + z). (5) 

A basic difference between the first and the second parameterization is that the later includes 
a possible epoch dependent correction which avoids the divergence at extremely high z. At 
low redshifts, when second order terms are neglected, the second parametrization reduces 
to the first one, that is, 77(2;) ~ 1 + r]QZ. Such one-parametric formulas arc also interesting 
because in the limit of extremaly low redshifts {z « 1), one finds ri{z) = 1 as should be 
expected since Dl — Da at this limit (see also parametrization II). In addition, for a given 
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data set, the likelihood of 770 must be peaked at r^o = 0, in order to satisfy the standard 
duality relation. In this context, by taking the SNe la from Constitution data (Hicken et al. 
2009) and galaxy clusters samples compiled by De Fillipis et al. (2005) and Bonamente et 
al. (2006), a direct test of the CDDR was acomplished (Holanda et al. 2010). As an extra 
bonus, the consistency between the strict validity of CDDR and the assumptions about the 
geometry based on elliptical and spherical /3 models was detailedly discussed. The sphericity 
assumption for the cluster geometry resulted in a larger incompatibility with the validity of 
the duality relation in comparison with an isothermal non spherical cluster geometry. More 
recently, such expressions were adopted by Gongalvez et al. (2011) in their studies on the 
validity of CDDR by including data from X-ray gas mass fraction of galaxy clusters. 

In this connection, it is also worth mentioning that Li, Wu & Yu (2011) rediscussed 
this independent cosmological test by using the latest Union2 SNe la data and the angular 
diameter distances from galaxy clusters thereby obtaining a more serious violation of the 
standard duality expression. In a simultaneous but independent work, Nair, Jhingan & 
Jain (2011) also investigated the strict validity of the CDDR by using the latest Union2 
SNe la data and the angular diameter distances from galaxy clusters, FRIIb radio galaxies 
and mock data. As an attempt to determine a possible redshift variation of the CDDR 
relation, they proposed six different (one and two indexes) parametrizations including, as 
particular cases, the ones adopted by Holanda et al. (2010, 2011a). As physically expected, 
their results depend both on the specific parametrization and the considered data sample. In 
particular, they conclude that the one index parametrizations, namely: rjv — rjg/il + z) and 
Vv — V9 exp{[2;/ (l + z)]/ (l + z)}, does not support the CDDR relation for the given data set. 
Meng, Zhang & Zhan (2011) also reinvestigated the model independent test by comparing 
two different methods and several parametrizations (one and two indexes) for ri{z). Their 
basic conclusion is that the triaxial ellipsoidal model is suggested by the model independent 
test at la while the spherical /3 model can only be accommodated at 3(T confidence level 
thereby agreeing with the results earlier derived by Holanda et al. (2010). 

It should be stressed that all the above described attempts to test the CDDR have been 
carried out based on a phenomenological approach. Usually, it is also not clear whether the 
nonstandard relation is the result of a modified luminosity distance or whether it should be 
associated to an extended angular diameter distance or both (see, however. Basset and Kunz 
2004). At this point, one may also ask whether such expressions to r]{z) can be derived from 
a more fundamental approach. In the affirmative case, it is also important to study their 
consequences for the present accelerating stage of the Universe. 

In this letter we consider both questions. Firstly, we show how any deformed CDDR can 
be derived by analysing possible theoretical modifications on the luminosity distance without 
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refraction effects. Analytical expressions for the dimensionless cosmic absorption parameter 
describing the above four parametrizations will be explicitly obtained in the framework of 
Gordon's optical metric as developed by Chen and Kantowski (2009a, 2009b) to include 
cosmic absorption. It will be also explicitly assumed that the angular diameter distances are 
not modified because their measurements involve only standard rulers and angular scales, 
and, more important, possible refractive effects have been neglected. Apart such hypotheses, 
the approach discussed here is quite general and can be applied for any deformation function, 
r]{z). Secondly, we also apply our results for the latest SNe la data. As we shall see, the 
modified luminosity distance can accomodate the observed supernova dimming even for a 
non-relativistic cold dark matter (CDM) Einstein - de Sitter model {Qm — 1)- The vahdity 
of the ACDM description is obtained only in the extreme limit of perfect cosmic transparency 
(negligible cosmic absorption). 



The concept of an optical metric was introduced long ago in a seminal paper by Gordon 
(1923). He proved the existence of a mapping between any solution of the general relativistic 
Maxwell's equations for a fluid with refraction index n{x) and the vacuum solutions in 
a related optical spacetime. In Gordon's treatment, only the refraction phenomenon was 
considered. More recently, by describing the Maxwell field as a monochromatic wave, Chen 
and Kantowski (2009a, 2009b) generalized such treatment in order to include the possibility 
of an absorption phenomenom in the Universe. As it appears, the presence of such an effect 
breaks naturally the validity of the standard distance duality relation as given by Eq. (1) 
since the light absorption violates the photon number conservation law. In this case, they 
prove that the luminosity distance takes the following form: 



where r denotes the optical depth associated to the cosmic absorption of the Universe and 
the superscript S specify the standard luminosity distance (no absorption) for which the 
Universe is assumed to be transparent. It is worth noticing that effects coming from a 
possible new physics like the interaction between photons and dark matter as discussed by 
Basset and Kunz (2004) are assumed to be neghgible (no new physics takes place). For a 
a spatially homogeneous and nondispersive (grey) absorption, the quantity r as derived by 
Chen & Kantowski (2009a) reads 



2. Luminosity Distance and Duality Relation 
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where a* is the dimensionless cosmic absorption parameter {a/Ho in the notation of Chen & 
Kantowski). In the above expression the refractive index was fixed to unity (n(z)=l), and, 
therefore, any possible refraction effect has been neglected. This is an important point, since 
in this case the standard angular diameter distances are not modified. 

Let us now discuss how the several functions r]{z) defining the deformed CDDR intro- 
duced in an ad-hoc way can be related with the dimensionless cosmic absorption parameter, 
a^:{z). As one may check, by inserting expression (6) into (2), we obtain that the cosmic opti- 
cal depth and the deformation function must be related by the simple expression e^/^ = r](z). 
In addition, this also means that any smooth r]{z) deformation function defines an effective 
cosmic absorption parameter given by: 



a.(z) = -^^ , (8) 

where a prime denotes derivative with respect to the reshift. Such a relation uniting the 
dimensionless cosmic absorption parameter and the deformation function, usually introduced 
by hand in the distance duality relation, is one of the main results of this section. When rj{z) 
is constant the cosmic absorption a^,{z) is identically null, and, therefore, r = 0. As remarked 
earlier, this also means that only the standard relation as determined by Etherington (1933) 
is possible at this limit, that is, rj^z) = 1 (note that r and rj are related by an exponential 
function). As one may check, the four deformed distance duality relations can analytically 
be expressed in terms of the dimensionless absorption parameter as: 

I. Basset and Kunz (2004): a^z) = 2[(/3 - l)^(^) + 7(1 + zy-% 

II. Avgoustidis et al. (2009): a^{z) = 2eE{z), 

III. Holanda et al. (2010): aJz) = + z)E{z) ^ 

1 + rjoz 

IV. Holanda et al. (2010): a^z^ - ^^^^^^^ 



1 + Z + TJoZ 

We see that the dimensionless Hubble parameter, E{z), and the free parameters appearing in 
the deformation functions define completely the cosmic absorption parameter as introduced 
by Chen and Kantowski (2009a, 2009b). 

At this point some comments are in order. Initially, we notice that in the limiting case 
7 = and (3 — e+1, the first expression for the cosmic absorption reduces to the second one. 
This should be expected since at such a limit the general deformation function of Basset 
and Kunz (2004) reduces to the one proposed by Avgoustidis (2009). Note also that such 
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expressions are satisfied for all Priedmann-Robertson- Walker (FRW) geometries (arbitrary 
values of Qk) and energetic content (baryons, dark matter, and dark energy). Since we are 
describing absorption, (a* > 0) this means that the parameters e and 770 must be positive. 

3. Extended Luminosity Distance and Supernova Dimming 

Nowadays, constraints based on SNe la data are considered the best method for studing 
the cosmic expansion history at z < 1.5. Let us now confront the extended luminosity 
distance including absorption as an optical metric phenomenon with the latest Supernova 
data. In our subsequent analyzes we consider only the fiat ACDM model for which the 
dimensionless Hubble parameter, E{z), takes the following form: 

E{z) = y/flM{l + z)^ + nA, (9) 

where — 1 — flM- 

To begin with, let us consider the Union2 supernova sample which is formed by 557 
measurements of distance moduli from Sne la as compiled by AmanuUah et al. (2010). In 
order to avoid effects from Hubble bubble, only 506 supernovae with redshifts greater than 
cz = 7000A;m/s were selected (Conley et al. 2007, Kessler et al. 2009, Sinclair et al. 2010). 
As widely known, the SNIa Union2 data are obtained by adding new datapoints (including 
the high redshift SNIa) to the original Union data (Kowalski et al. 2008). For this enlarged 
sample, a number of refinements to the original Union analysis chain has been done, in 
particular, the relative importance of systematic effects was higlighted (in this connection 
see also Sullivan et al. (2011) to the Legacy Survey sample (SNLS3)). 

In our statistical analysis we consider a maximum likelihood determined by a statis- 
tics 

SNIa ^obs,i 

where Hthi^i) = 5 logj^g ^Li^i, p) + iiq is the theoretical distance moduli, /Iq — 25 — log^g -^O; 

is the luminosity distance from Eq. (6), aobs,i is the uncertainty in the individual 
distances (including systematic errors), and the complete set of parameters is given by 
p = {Ho,Qm:(x*)- It should be stressed that for the considered SNe la subsample we 
have combined the statistical plus systematic errors in quadrature as compiled in Table 7 
of AmanuUah et al. (2010) and neglected af^. We have also marginalized on the Hubble 
constant. 

In Fig. 1(a) we show the contours on the Qm — plane corresponding to a flat ACDM 
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model and by considering that the absorption a* is constant. As indicated, the shadow hnes 
are cuts in the regions of 68.3%, 95.4% and 99.7% of probability. The constraints on the 
free parameters are restricted to 0.0 < a* < 1.55 and 0.21 < Qm < 1-0 at 2a of statistical 
confidence while to the latter we have obtained 0.0 < a* < 1.46 and 0.22 < Qm ^ 1-0 at 
2a of statistical confidence. In the absence of absorption {a^, = 0) the limit on the density 
parameter is 0.21 < Qm < 0.34 (concordance model) while for Qm = 1 (Einstein de Sitter), 
we find that the absorption parameter lies on the interval 1.20 < a* < 1.54 {2a). Note that 
the positiveness of a* (absorption), implies that a pure de Sitter model — 1, Qm — 0) is 
not allowed by such data. The best fit to the free parameters are Qm — 0.27, a* = 0.0 with 
xLn = 330.5. 

In Fig. 1(b) we display the likelihood distribution functions, e~^^/^ of a*, for the 
Einstein-de Sitter model {Qm — 1) with constant absorption. The blue curve includes 
statistical + systematic errors, but, for comparison, we have also shown the red curve with 
only statistical errors. The upper and lower horizontal lines correspond to 68.3%, and 95.4% 
C.I., respectively. By marginalizing on the Hubble parameter we obtain a* = 1.38±0.08(0.15) 
with la (2a) of probability and Xmin — 331.3. In this model the Universe is always decelerat- 
ing since q{z) = qo = 1/2. Instead of dark energy we have a cosmic medium whose absorbing 
properties is quantified by the dimensionless parameter a* ~ 1.4 which is responsible for the 
SNe la dimming. The dimensional absorption parameter, a — a^^Ho ~ 10~^Mpc~^, is nearly 
the same one previously obtained by Chen and Kantowski (2009a). 

Let us now consider the deformation function, rj^z) = (l+z)"^, as adopted by Avgoustidis 
et al. (2009). As shown in the previous section, the associated cosmic absorption parameter 
in this case is a* (2;) = 2eE{z). In order to simplify the notation, in what follows we consider 
the parameter, ccq = 2e. 

In Figure 2(a) we display the corresponding plots for the {Qm,(^o) plane to the case 
of variable absorption, a* = aoE{z). As in Fig. 1(a) the analysis includes statistical plus 
systematic errors. The confidence region (2a) in this plane is defined by 0.0 < ckq < 0.92 
and 0.20 < < 1-0. Again, the best fit is the ACDM with best fit = 0.27, ao = 0.0 
with xLn = 330.5. 

In Fig. 2(b) we display the likelihood distribution functions for ckq. Blue and red curves 
correspond to the same analysis of Fig. 1(b), that is, with and without systematics. To the 
latter case, we obtain that ckq = 0.88 ± 0.05(0.10) with la {2a) of confidence level. Once 
again we have q{z) = go = 1/2 in the presence of the absorption, but with Xmin ~ 335.9. 

It should be stressed that the above results for the Qm — and VLm — cto planes, based 
on a subsample of the Union2 SNe la data, are strongly suggesting that one of the following 
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possibilities must be true: 

(i) The Universe has an extra dark energy component (A) which provides the fuel for 
the observed accelerating stage (go < 0,a^ = 0). 

(ii) We live in a decelerating Einstein-de Sitter Universe {Qm = 1) endowed with a 
cosmic absorption mechanism which is responsible for the dimming of the distant Supernovae 
la (go = 1/2, a, ^0). 

(iii) We are not living in the extreme cases above mentioned, that is, both phenomena 
(cosmic oppacity and A) are partially operating in the observed Universe. 

In conclusion, by adopting the Gordon-Chen-Kantoswski description we have shown how 
to obtain analytical expressions for any deformed distance duality relation (parametrized by 
r){z) function) in terms of the cosmic absorption parameter and vice- versa. Four especific 
examples of r]{z) functions, recently proposed in the literature, were explicitly considered. 
The associated absorption parameter for each one was explicitly determined and some of 
its physical implications also were discussed in detail. Two different scenarios with cosmic 
absorption were proposed and their free parameters constrained by using the SNe la data 
(AmanuUah et al. 2010). Working in the in the framework of a flat ACDM model we 
have found that both scenarios (a* = aoE{z) and a* = constant) are also compatible with 
a decelerating Einstein-de Sitter Universe (see Figures la and 2a). Finally, as remarked 
before, in the present study we have neglected any effect due to a possible light refraction 
phenomenon as described in the original Gordon's description. The general case it will be 
discussed in a forthcoming communication. 
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Fig. 1. — a) Contours in the Qm — C(* plane for 506 supernovae data by considering a flat 
ACDM model with constant absorption parameter, a*. Note that the Einstein-de Sitter 
model is now allowed by such data, b) The likelihood function of the constant absorption 
a^, for the Einstein-de Sitter model {Qm = !)• The red line takes into account only statis- 
tical while the blue curve includes both errors (statistical + systematic). The upper and 
lower horizontal lines correspond to 68.3% (Icr), and 95.4% {2a) c.l., respectively. The free 
parameter is constrained to a* = 1.38 ± 0.08 (statistical + systematic, la c.l.). 
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Fig. 2. — a) Contours in the Qm ~(^o plane for the same 506 supernovae data by considering 
a flat ACDM model with an epoch-dependent absorption, a* = aQE{z), as fixed by the 
Avgoustidis et al. (2009) deformation function (see sections I and II). The best fit to the 
pair of free parameters is given by {Qm, ao) = (0.27, 0.0). b) The likelihood functions of ao in 
the case of fiat Einstein-de Sitter model. The free parameter is constrained to = 0.88±0.05 
(statistical + systematic, la c.l.). 



